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■ The decays — > 27r+27r and Itt^tt^tt are studied using various 



effective Lagrangians for tt and p (and in some cases ai) mesons, all of which 
respect the approximate chiral symmetry of the strong interaction. Partial 
. widths of the order of 1 keV or less are found in all cases. These are an order 

in 

. of magnitude smaller than recent predictions based on non-chiral models. 

o 

I There have been recent hopes that the rare decays —>■ 27r+27r^ and 27r*^7r+7r~ 

might be observable in the near future in experiments at high-luminosity e~^e~ machines, 
^ I such as VEPP-2M ||^ or DA$NE [Q. These decays might test various aspects of the low- 
energy effective Lagrangians that have been proposed for the interactions of p-mesons and 
^ I pions [§-01 • For example, the 27r°7r+7r mode is sensitive to the 3p coupling that appears 
in theories with gauge-type couplings of the p. Current experimental limits on the partial 
widths for the decays are 30 keV for the 27r+27r^ mode ^ and 6 keV for the 27r'^7r+7r~ mode 
0. These are already sufficient to rule out some earlier estimates which did not incorporate 
chiral symmetry, for example one dominated by vroi and 7ra2 intermediate states [l^ . 



Bramon, Grau and Pancheri studied the 27r+27r~ decay with two commonly-used 
formalisms for including the p meson in low-energy effective Lagrangians. For the simplest 
"hidden-gauge" theory of Bando et al. [pT| , p!2| they obtained a partial width of 7.5 ±0.8 keV. 



In contrast, for a simplified "massive Yang-Mills" theory, where the p is coupled to a sigma 
model as an SU(2)y gauge boson they obtained 60 ± 7 keV. This suggested that the 



process could distinguish between these approaches, and indeed that the massive Yang-Mills 
one is already ruled out by the data. 

More recently, Eidelman, Silagadze and Kuraev have pointed out that the massive 
Yang-Mills model of Ref. does not respect chiral symmetry. One possibility is to use 
gauge couplings corresponding to the full SU(2)j:jxSU(2)j;^ symmetry by introducing the ai 
as the chiral partner of the p [ll3|,|5|. An alternative, and the one followed in Ref. P|, is 



that proposed by Brihaye, Pak and Rossi where additional terms are added to the naive 



TT-p Lagrangian in order to preserve low-energy theorems. Including Air and pAn correction 
terms, the authors of Ref. obtained partial widths of 16 ± 1 keV for the 27r+27r~ mode 
and 6.0 ± 0.2 keV for the 27r'^7r+7r~ mode, consistent with the present experimental limits 
and within reach of future experiments. 

Although the Lagrangian of Ref. is constructed to satisfy vector-meson dominance 
and a range of low-energy theorems for vrvr scattering and the couplings of photons, it is still 
not chirally symmetric. A Yang-Mills Lagrangian for pions and p mesons only contains a 
momentum-independent 27r2p coupling of the form 

^gUiL X p/. (1) 

This leads to a vr-p scattering amplitude that does not vanish at threshold in the chiral limit. 



and so is in conflict with low-energy theorems of chiral symmetry |T^. One can add further 
terms to such a Lagrangian to cancel out this unwanted behaviour, as described in Ref. [0. 
It is, however, straightforward to work with a manifestly chiral approach from the start. We 
have therefore calculated the p^ — > Arc decay rates using several such Lagrangians. 

The low-energy interactions among pions can be described by a chirally symmetric sigma 
model. There is, however, considerable freedom in extending such a description to incor- 
porate spin-1 mesons 1^-0]. We have examined three such approaches: i) the hidden-gauge 



theories of Bando et al. [0,0 in which a local SU(2) symmetry is introduced into the non- 
linear sigma model and the p is treated as the gauge boson of this symmetry; ii) the massive 
Yang-Mills theories in which p and ai fields transform linearly under the chiral symmetry; 



iii) the approach suggested by Weinberg and developed by Callan, Coleman, Wess and 



Zumino (CCWZ) in which the spin-1 mesons transform homogeneously under a non- 
linear realisation of chiral symmetry. One should remember that, in general, all of these 
are equivalent and any Lagrangian of one approach can be converted into a corresponding 
Lagrangian of another 

If we work at tree level, the relevant diagrams are shown in Figure 1. Including anomalous 
terms in the effective Lagrangian leads to additional contributions, similar to those of 1(f) 
and 1(g) but involving an intermediate lj meson. In all cases we include a symmetry-breaking 
term of the form 

&mlTriU + U^), (2) 

where U = exp(zr.7r//^). This term produces the non-zero pion mass. It also has the effect 
of modifying the 4tt vertex in diagram 1(b). This is quantitatively important for our results, 
since they are much smaller that those of Refs. |jl|,§] , but it does not affect their qualitative 
features. 

First, we consider the simplest example of a hidden-gauge theory, defined by the La- 
grangian of Bando et al. [|ll|]. To allow for contributions from the anomalous part of the 
Lagrangian, we introduce the u) meson by enlarging the gauge group to f/(2)y. There are 
then three possible gauge- covariant anomalous terms, with undetermined coefficients. We 



adopt the suggestion of Ref. |18| with regard to these coefficients, including an upn vertex 



but no wStt contact term. The parameters of this model satisfy the relations 

ml = a-a'^fl, QpTTn = ^g- (3) 



The choice a = 2 yields the KSRF relation |T9[, universal coupling of the p and allows 
vector dominance to be realised in the 77r7r coupling. In our calculations we take the following 
values for the parameters: 

= 92.4MeV, = 139.6MeV, = = 770MeV, a = 2, g = 5.89. (4) 

The amplitudes for the p° — > Att decays are straightforward to derive, but the final 
expressions are rather long and so we do not present them here. Having found them, we 



must integrate over phase space to obtain the corresponding partial widths. We express 
the integrals in terms of the variables of Ref. I^T] and evaluate them numerically using the 
NAG routine DOIFDF, which maps the region of integration onto an ra- dimensional sphere 



and uses the method of Sag and Szekeres |22| to perform the integration. With a suitable 
choice of the parameter controlling the mapping onto the sphere, we find that about 50,000 
integration points are sufficient to give the integrals to an accuracy of one part in a thousand. 

The results obtained in this model are shown in the first line of Table 1, labelled HG. 
They are roughly an order of magnitude smaller than any of the results given in Refs. 
Bramon et al. [Q have calculated the process p° 27r+27r^ in the same model and they quote 
a value of 7.5 ± O.SkeV, in sharp contrast with our result. The crucial difference between 
the calculations lies in the strength of the direct p47r coupling. Bramon et al. assumed that 
the expression for this vertex was the same as in a massive Yang-Mills model, specifically 
the p47r term of 



In fact, the corresponding term in the hidden-gauge model should be written, in the unitary 
gauge, as 

-2tagf^T:t[p^{u^d^u + ud^u^)] = ^9 {^^ - ^2l' + ■ ■ P^-K x 9^ir. (6) 

where u is the square root of U. Although both expressions yield the same pmi coupling, 
the p47r terms differ by a factor of four. Hence one cannot take the latter coupling to be 
the same as in a massive Yang-Mills model. Reducing the contribution of diagram 1(a) by 
a factor of four has a large effect on the total amplitude, explaining the difference between 
the result of Ref. [|I| and ours.Q 



^Comparison of the amplitudes given by Refs. [T|,^ also indicates a discrepancy in the evaluation 
of diagram 1(b). Our evaluation concurs with that of Ref. [0. However, the impact of this error is 
fairly small. 
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We have examined the importance of the anomalous processes for our results. Keeping 
only the non-anomalous contributions leads to the result for —>■ 27r'^7r+7r^ labelled HGNA, 
showing that the two types of term are of similar importance. We have also looked at the 
effect of omitting the symmetry-breaking An interaction of Eq. (2) (but still keeping the 
physical pion mass in the propagator etc.). The results, labelled HGCS, indicate that this 
term is a smaller, but still significant, contribution. As a simple estimate of contributions 
beyond tree level, we have examined the effect of including the finite width of the p in its 
propagator (as in Ref. f^). We find that this does not alter the results significantly. 

We now turn to the massive Yang-Mills type of theory, so named because the spin-1 
mesons are introduced as if they were gauge bosons of chiral symmetry. Local chiral sym- 
metry is broken by the mass terms for those mesons. As mentioned above, chiral symmetry 
requires that the ai be included as the chiral partner of the p in these models. We start 
with the simplest example of such a Lagrangian, which is just a version of that in Refs. 
23| , |2^ but using a non-linear rather than a linear sigma model. In this minimal model the 



p and Oi masses satisfy Weinberg's relation, m^^ = \/2mp |2^. Note that, in this approach, 
chiral symmetry requires strong cancellations among the contributions to any amplitude. 
One cannot replace by hand, say, the a\ mass by its physical value in such a calculation 
without violating chiral symmetry. The constraints imposed by the symmetry mean also 
that one must not omit the a\ field from this approach without introducing compensating 
terms that maintain chiral low-energy theorems 0,0] ■ 

A feature of this approach is the appearance of a vr-ai mixing term in the Lagrangian. To 
remove this and diagonalise the free-field part of the Lagrangian, one can subtract a term 
proportional to 9^7r from the axial field and renormalise the pion field. This constitutes 
a minimal diagonalization procedure and leads to physical fields with complicated chiral- 
transformation properties, but it is sufficient for our present purposes. The procedure gen- 
erates various additional three- and four-point interactions arising from the gauge-invariant 
kinetic energies of the spin-1 fields. These include a pTrvr interaction of order Oijf'^ which 
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reduces the pun coupling strength at mp by a factor of | compared to its value at zero 
momentum. The minimal model is thus unable to give a good description of the p and 
Oi masses and widths. To remove this deficiency, the model must be supplemented with 
additional terms |T3|,|5|, a point we return to below. 

For the moment though, we consider the minimal model with the physical value of rrip 
and g = Qpt^t, = 5.89. We extend the model to include the u meson, taking the relevant 
anomalous upn and cuSvr vertices from Ref. The results for p° — » An are shown in 

Table 1, labelled MMYM. They are similar in magnitude to those of the hidden-gauge 
model. As already mentioned, the calculations of Refs. used Yang-Mills models that do 
not respect chiral symmetry. Without the ensuing cancellations, they lead to partial widths 
that are an order of magnitude larger than ours. 

Finally we consider several Lagrangians of the CCWZ type |1T7|J6||. These are expressed 
in terms of spin-1 fields that transform homogeneously under a non-linear realisation of 
chiral symmetry [|1^,|T^. Our reasons for this are twofold. First, the models just described 
can be converted by a change of variables into equivalent CCWZ Lagrangians 0,0, which 
should yield the same predictions for any observable as the original models. These therefore 
provide a useful check on our results above. Second, the CCWZ formalism provides a 
convenient framework to examine the sensitivity of our results to assumptions about the oi. 
In contrast to the massive Yang-Mills approach, the parameters describing the ai mass and 
couplings may be independently changed without the need to introduce further terms into 
the Lagrangian. 

Apart from labelling of the coefficients, we use the notation of Ref. to express the 
relevant terms in the non-anomalous Lagrangian as: 



6 



+ lc3Tr ((K,K] - K,K])') - ^C4Tr(K,«J(K,An - [u\A^'])). (7) 

The minimal Yang-Mills model discussed above corresponds to the following set of coeffi- 
cients: 

9i = ^{'^-Z^), 92 = '^9, 5-3 = 5-4 = ^^ 

Cl =5-1, C2 = 1 - C3 = Z^, C4 = 9i93, (8) 



where 



2 -f 2 

Z^ = l- (9) 



The appropriate coefficients for the simplest hidden-gauge model |jTT| can be obtained by 
setting Z = in Eq. (8) and dropping all terms involving the ai. We have calculated the 
amplitudes for p° 4tt within this framework and verified the results for the non-anomalous 
parts of the hidden-gauge and massive Yang-Mills models described above. 

The anomalous sectors of these models can also be converted into CCWZ form, although 
this is somewhat involved for the massive Yang-Mills approach. However the sum of am- 
plitudes for the anomalous diagrams must remain unaltered by the change of variables. For 
simplicity, therefore, we take the anomalous part of the amplitude directly from the original 
version of the massive Yang-Mills model. For the couplings given in Eqs. (8,9), rrii^ = 783 
MeV and rria^ = 1230 MeV, we get the results labelled MYM+1 in Table 1. These corre- 
spond to a massive Yang-Mills Lagrangian with non-minimal terms of the type suggested in 
Refs. igjl. 

The parameter choice of Eqs. (8,9) still suffers from the fact that it gives too small a 
width for p 27r. It is straightforward to modify the CCWZ parameters to remove this 
deficiency. This is equivalent to adding a further non-minimal term to the massive Yang- 
Mills Lagrangian ||l3i|5| in order to cancel the diagonalization-induced 0{p^) pnir coupling. 
In the CCWZ representation, the corresponding parameters are 

9i = 7^, 92 = 29, 93 = 0, 5(4 = ^^ 
2^ 

Cl = -i^(l-^'), C2 = l, C3 = Z\ C4 = 0. (10) 
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The results for this set, with the empirical meson masses, are labelled MYM+2 in Table 1. 

We have examined the sensitivity of our results to the 3p coupling present in these 
Lagrangians. In the hidden-gauge and Yang-Mills models described above this coupling is 
equal to the ptttt one because of the assumed universal coupling of the p. Using the CCWZ 
equivalents of these models, we have varied g2 by ±30% and found that the decay rate for 
p° 27r'^7r+7r^ changes by about ±1%. 

These results show that, for all of the chirally symmetric models considered, the partial 
widths for the decays 27r+27r^ and 27r°7r+7r^ are of the order of 1 keV, corre- 

sponding to cross sections of the order of 5 pb. These are an order of magnitude smaller 
than the predictions of Refs. Although the processes may be hard to observe in fu- 

ture experiments, they should not be beyond the reach of DA$NE which is designed to 
have a luminosity of 5 x 10*^ b^^s^^ The precise results are sensitive to the choice of 
Lagrangian and parameters, receiving significant contributions from anomalous processes 



(about which there is still some debate |^) and symmetry-breaking interactions. However, 
any uncertainties in the determination of the relevant parameters are not sufficient to affect 
our general conclusion on the magnitudes of the decay rates for for p° ^ 47r . 
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TABLE 



Model 


pO _^ 27r+27r- 


pO _^ 27r°7r+7r- 


HG 


0.89 


0.44 


HGNA 


0.89 


0.24 


HGCS 


0.59 


0.37 


MMYM 


0.68 


0.37 


MYM+1 


0.63 


0.34 


MYM+2 


1.03 


0.39 



Table 1. Decay widths (in keV) for p° — > Air in various versions of the chiral models. See 
the text for the definitions of these. 
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FIGURE 




Figure 1. Diagrams contributing to the p° An decays in chiral effective Lagrangians of tt, 
p and oi mesons. Single lines denote pions, double lines spin-1 mesons. Anomalous terms 
in the effective Lagrangian introduce diagrams similar to (f) and (g) but with an cu meson 
replacing the ai. 
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